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We performed a finite-temperature quantum Monte 
Carlo simulation of the one-dimensional spin-| Heisenberg 
model with nearest-neighbor interaction coupled to Einstein 
phonons. Our method allows to treat easily up to 100 phonons 
per site and the results presented are practically free from 
truncation errors. We studied in detail the magnetic suscepti- 
bility, the specific heat, the phonon occupation, the dimeriza- 
tion, and the spin-correlation function for various spin-phonon 
couplings and phonon frequencies. In particular we give ev- 
idence for the transition from a gapless to a massive phase 
by studying the finite-size behavior of the susceptibility. We 
also show that the dimerization is proportional to g 2 /fi for 
T <2J. 



PACS numbers: 75.10 Jm, 75.40.Mg, 75.50.Ee 



I. INTRODUCTION 



In this paper we study the finite-temperature proper- 
ties (T/J > 0.1) of the spin-i Heisenberg model coupled 
to (dispersionless) Einstein phonons, 
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JV 



i=i 

N 

+fi$>fa, (i) 

i=i 

in the parameter range of the spin-phonon couplings < 
g < 1.5 and phonon frequencies 0.2 < fi/J < 20. 

The phenomenological interest in .-this model dates 
back to the fundamental work of PytteEl who showed that 
a spin-i Heisenberg chain coupled to three-dimensional 
phonons undergoes a dimerization transition, generally 
referred to as. "spin-Peierls" transition. Motivated by 
the discoveryB of a transition of this type in the inor- 
ganic substance CuGe03 , the model proposed by Pytte 
and its modifications have been widely studied in recent 
time. (Thus following a period when CuGeOa was mostly 
discussed in the framework of a frustrated and dimerized 
spin-Hamiltonian with a temperature-dependent dimer- 
ization as a relict of the spin-phonon coupling.) In partic- 
ular the low-temperature properties of the model Eq. (pi) 



were evaluated recently by Lanczos and DMRG tech- 
niques in Refs. and by a continuous time world-line 
Monte Carlo algorithm in Ref. ||. At T ^ a stochastic 
series expansion Monte Carlo was performed in Ref. [y. 

In this paper we are mostly concerned with the finite- 
temperature properties of the model. We thus partly 
reproduce results of Ref. ^| where the model Eq. (|l|) was 
studied at finite T for g = l/(4\/2) and = 0.2 J. In 
contrast to Ref. ^ we analyze in this work both the adia- 
batic limit of the model and also its gapless region where 
the thermodynamic observables are only slightly effected 
by the coupling to the phonons. We also calculate the 
specific heat, the spin-spin correlations and take advan- 
tage of the pronounced periodic patterns in the local 
phonon displacement to determine the parameter range 
where the model is dimerized. Our choice of parameters 
is motivated by Ref. [l(] where the dominating phonon 
modes for CuGeOa are found to be fli/J = 2, H2/J = 4, 
g x = -1/8, and g 2 = 1/2. 

The outline of the paper is as follows. In Sec. [n] we 
briefly describe the Monte Carlo loop algorithm we em- 
ployed. In Sec. [II we present the thermodynamic proper- 
ties, in particular we demonstrate the presence of a phase 
transition from a gapless to a gapped phase by studying 
the finite-size behavior of the magnetic susceptibility. We 
also show the effect of the spin-phonon coupling on the 
specific heat and discuss results for the phonon occupa- 
tion number. In Sec. |y| we consider the local phonon 
displacement operator and determine the approximate 
boundary of the dimerized phase of the model. We study 
the dependence of the effective coupling J c g on fl and 
g and finally find an interesting crossover in the spin- 
spin correlations between the high and low-temperature 
phase. Conclusions are given in Sec. 0. 



II. QUANTUM MONTE CARLO METHOD 

For this work we developed an extension of the quan- 
tum Monte Carlo loop algorithmlidO. Like for the world- 
line method the key idea of the loop algorithm is a map- 
ping of the partition function of a n-dimensional quantum 
spin system onto a (n+l)-dimensional classical system. 
However, its substantial advantage is that compared to 
the world-line method it allows global spin updates, thus 
substantially reducing the autocorrelation times. Also 
the variation of both winding number and magnetization 
are automatically included. 
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For our simulation we used two Trotter layers for the 
spin-spin and the spin-phonon part of H and introduced 
a third layer for the free phonon part. We then com- 
bined global loop updates for the spin-phonon and spin- 
spin coupling with local heat bath updates modifying the 
phonon occupation number. It is obvious that detailed 
balance is satisfied for the two steps and thus also for 
the whole procedure. The decomposition of the partition 
function which forms the basis of our calculation is given 
in the appendix. 

We studied chains with 16 up to 128 sites using typi- 
cally 10 5 . . . 10 6 spin updates for each temperature. Even 
though the importance sampling technique was employed 
for the phonon update, its acceptance rate is still smaller 
than that of the spin update. Best results were obtained 
by making « 20 phonon updates per spin update and us- 
ing only the last of the phonon updates for the evaluation 
of the expectation values. Also for each temperature the 
first 10% of all the sweeps were used for thermalization. 
We find that the measured quantities depend linearly on 
the inverse Trotter number squared. This can be seen in 
Fig. 0, where we show the susceptibility for a N = 16 
system and Trotter number M = 10, . . . , 80 as function 
of 1/M 2 . In the present work we either extrapolated our 
results in the Trotter number, according to the above law, 
or when we found the fluctuations to be larger than the 
effect of the finite Trotter number, we gave the explicit 
value of the Trotter number in the figures. 
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FIG. 1. Magnetic susceptibility for temperatures 

T/J = 0.2, 0.5, 0.8, 1.1, 1.4 as a function l/(Trotter number) 2 
for N — 16 spins, phonon frequency £1 = J, and spin-phonon 
coupling constant g = 0.2. 



III. THERMODYNAMIC PROPERTIES 

A. Magnetic Susceptibility 

It is known that the magnetic susceptibility of CuGe03 
in the temperature range above the spin-Pcierls tran- 



sition j-eaa, be well fitted by a frustrated Heisenberg 
mode£3~lij. It is often objected, however, that the ex- 
cellent agreement of theory and experiment is acciden- 
tal, since it does not take into account interchain inter- 
actions and spin-phonon couplings. Motivated by this 
controversy we pursue the question to what extent the 
spin-phonon coupling influences the susceptibility. 

We summarize our observations on the high- 
temperature properties of the susceptibility in the fol- 
lowing three points. 

(i) For fixed phonon frequency, the overall height of the 
susceptibility is lowered with an increasing spin-phonon 
coupling (Fig. |). 

(ii) For fixed g, the susceptibility is growing and its 
maximum is shifted to lower temperatures with increas- 
ing frequency fl (Fig. ||). 

(iii) As the spin-phonon coupling is increased the max- 
ima of the susceptibility curves are shifted to higher tem- 
peratures. 
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FIG. 2. Magnetic susceptibility versus temperature for 
£1 = 2 J and spin-phonon coupling g between 0.2 and 1.5. 
The number of sites is N = 64 for low temperatures, T < J, 
and N = 16 for high temperatures. For comparison we show 
the exact-results for the Heisenberg model calculated by A. 
KlumperES via the quantum transfer matrix approach. 

Point (i) is not surprising, since an increase of the spin- 
phonon coupling is expected to decrease the magnetic 
order. Point (ii) can be explained by the fact that for 

3> gJ the magnetic and phononic degrees of freedom 
decouple. So for Q — > oo (and gJ fixed) the magnetic 
properties are again dominated by the antifcrromagnetic 
Heisenberg model. Point (iii) is of some phenomcnolog- 
ical consequence because it means that for a substance 
where the spin-phonon coupling cannot be neglected, the 
value of the spin-spin coupling is overestimated, if deter- 
mined by comparing the experimental data with the sus- 
ceptibility calculated in a pure spin model. (We assumed 
that J is determined by matching temperatures where 
the experimental and theoretical susceptibility curves 
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have their maxima, which is a standard procedure, see 
for example Ref. 

Also it should be noted that for frustrated models 
T m axl J is shifted to lower temperatures with increasing 
frustration parameter, i. e. the effect of the spin-phonon 
coupling is opposite to that of frustr ation. Together with 
the argument given in Section IV B it is obvious that the 
susceptibility of CuGeC>3 cannot be fitted with a model 
of the type Eq. (Eh. 



dimerization patterns at low T. 

The existence of this phase transition has been previ- 
ously reported at T = in Ref. || where by means of a 
continuous time world-line algorithm it is shown that a 
critical coupling between gapless and spin-fluid phase ex- 
ists for a fixed fl. Also for a model with a slightly different 
spin-phonon coupling the phase diagram was determined 
in Ref. E_7] by means of DMRG technique and for a model 
including a frustrating spin-spin coupling evidence for the 
transition can be found in Ref. [s]. 
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FIG. 3. Magnetic susceptibility versus temperature for 
g — 1.0 and phonon frequency between Q = J and Q = 20J. 
The number of sites is N = 64 for low temperatures, T < J 
(T < 1.5J for the fi = J chain), and N — 16 for high rtem- 
peratures. For j = 0we show the exact result from Ref.cJ. 

To draw definite conclusions from the susceptibilities 
at low T, finite-size effects have to be studied carefully. 

In Fig. |a we show results for £1 = 2J and g = 0.5 
and chain lengths N — 16,32,64, 128. The tendency of 
the susceptibility curves is to approach a finite value at 
T — > 0. Thus the sharp decrease in the curves for N = 
16, 32, 64 is a finite-size effect, i.e. the energy gap causing 
the drop is due to the discreteness of the spectrum of the 
finite system, it is not an intrinsic property of the spin 
chain (Fig. 

The situation is drastically changed if we consider 
g = 1.5 and O = J (lower curve in Fig. Here the 

drop is unchanged when going to larger systems. This 
means that the gap in the spectrum is larger than the 
gap generated by the finite size of the chain. The upper 
curve in Fig. ^b shows results for g — 0.8 and f2 = J, 
which exhibit clear finite-size effects but where the na- 
ture of the gap cannot be determined from systems of 
size N < 128. Since for g — the model is gapless and in 
the adiabatic limit the model is known to be massive, the 
finite-size behavior shown in Fig. ^b gives clear evidence 
that the transition between the two phases occurs at fi- 
nite g. This transition to a gapped phase is a signature 
of the lattice distortion. Indeed we will show in section 



IV A, that in the gapped phase we find well developed 
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FIG. 4. Finite-size behavior of the magnetic susceptibil- 
ity. The upper figure shows a chain with phonon frequency 
f2 = 2 J and spin-phonon coupling constant g = 0.5, the lower 
figure shows chains with Q = J and g = 0.8 (upper curve) 
and g — 1.5 (lower curve). 



B. Specific Heat, Phonon Occupation Number 

The idea to get further information about magnetic 
properties by considering besides the susceptibility also 
the specific heat suggests itself since the specific heat is 
an observable, that can be as easily measured and cal- 
culated as the susceptibility. However, in practice only 
for substances with vew small J (Cu(L-alanine)2 with a 
J = 0.55K for exampleOJ) the specific heat can be directly 
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compared with experiment. For most of the inorganic 
substances (in particular for CuGe03 ) which have spin- 
spin couplings of the order of 100K, lattice vibrations 
dominate the specific heat and it is an unsolved problem 
to extract its magnetic part. As a first approximatior£3 
it was assumed for CuGeC>3 that it is sufficient to sub- 
tract the contribution of the phonons and consider the 
rest as the contribution of the spins. We show that such 
an approximation is not justified even in the simple case 
of dispersionless phonons. 
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oretical results in a quantitative way with experimental 
data. 

Further insight in the model can be gained, by studying 
to what extent the phonon-occupation number is influ- 
enced by the spin-phonon coupling. 

1 N 

^) — E(^}- (3) 

i=i 

We first consider fi = 4J and g = 1, parameters for 
which we found above that the specific heat can be ap- 
proximated without the spin-phonon coupling. For n(T) 
we find in this case that the effect of the spin-phonon 
coupling is solely to shift the contribution 

n f^ = ^WZTJ ( 4 ) 

of the free phonons by a constant value no (see Fig. 
We determined from the low temperature data that no = 
2(gJ/$Y) 2 , which in agreement with Ref. where this 
relation was found for a chain with a single k = n phonon 
mode. However, for g = 1 and Q/J = 1,2 deviations 
from the free phonon case become visible. In particular 
the slope of n(T) at high T is substantially smaller than 
l/f2. Also the constant shift deviates from 2(<?J/f2) 2 . 
Note that the drop of n(T) at low T is an effect of both 
finite Trotter number and the finite size of the system. 



FIG. 5. Temperature-dependence of the specific heat for 
a chain with phonon frequency Q. — J and spin-phonon cou- 
pling constant g — 1.0. The number of sites is 16 and 64, 
respectively. 

To demonstrate this we plot in Fig. || the specific heat 
c for a chain with il = J and g = 1.0 together with the 
specific heat c/ of the free Einstein phonons 

/tt\ 2 e Q/T 
Cf = [t) (e^TTf' (2) 

with = J and— the contribution of the isotropic 
Heisenberg modefel As can be clearly seen the Einstein 
phonons are dominant at high temperatures, in particu- 
lar the constancy of the specific heat, as expected accord- 
ing to the Dulong-Petit rule, is visible. However, also the 
effect of the spin-phonon coupling can be perceived, since 
we find that the maximum at T J of the Cf(T)+cn(T) 
(dot-dashed line) becomes less pronounced and is shifted 
toward higher temperatures by the presence of the spin- 
phonon coupling. When going to larger frequencies (and 
g fixed), the contribution of the spin-phonon coupling de- 
creases, but only for > 4J and g = 1.0 we find that 
c(T) can be approximated by c/(T) + ch(T). Of course 
when considering the specific heat the phonon dispersion, 
which is not taken into account, is most important and 
more work on realistic models is needed to compare the- 
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FIG. 6. Mean phonon occupation number as a function of 
temperature for phonon frequency Q = 4J and spin-phonon 
coupling constant g = 1.0. The number of sites is 16 and 
64, respectively. The solid line shows n/(T) + hq(T), where 
no(T) = n(T) — rif(T) was approximated by the tempera- 
ture-independent value no(0) = 2(gJ/Q) 2 . 

We conclude this section with a more technical re- 
mark concerning truncation errors. As the dimension 
of the subspaces is quickly growing with the number of 
phonons the limits of a method requiring truncation is 
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soon reached, this is relevant when using exact diago- 
nalization which allows to treat only « 16 spins and a 
few phonons per site. At zero temperature the trun- 
cation problem can-be partly overcome by using coher- 
ent phonon statesa€2l or employing the DMRG technique, 
which allows to treat system sizes of up to 256 spinsEJ. 
At finite temperatures, however, the truncation problem 
is even more important. With the probability distribu- 
tion of the phonon occupation number shown in Fig. |?j 
we dramatically demonstrate that high cutoffs are essen- 
tial when studying finite temperature properties. For the 
correct treatment of a chain with Q > J and g = 1.0 at 
T = J a truncation at ~ 30 spins is needed (as can be 
seen from the figure the probability for more than 28 
phonons per site is less than 10~'). For 0,/J — 0.5 the 
truncation is < 50 spins/site and only very high phonon 
frequencies as il/J = 20 can be treated with 4 phonons 
per site. 
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FIG. 7. Probability of the phonon occupation number for 
phonon frequencies between Q./.J — 0.5 and 20 for a chain 
with N — 64 sites, Trotter number 80, spin-phonon coupling 
constant g = 1.0, and temperature T = J. 



IV. DIMERIZATION AND 
SPIN-CORRELATIONS 

Addressing now the question of the dimerization, we 
expect that the ground state of the model Eq. (||) in the 
massive phase shows short range correlations and strong 
signatures of dimerization. We extract the boundaries 
of the phase diagram in the g-£l plane by studying the 
dimerization patterns at low temperatures and corrobo- 
rate the results by considering the divergence of the struc- 
ture factor at q = w. We also give account of the averaged 
phonon displacement often referred to as an "effective 
coupling" , which is nonzero for all g ^ and Q ^ 0. 



A. Dimerization 

To begin, let us recall the definition of the local phonon 
displacement operator 



S^gibl + h) 



(5) 



which corresponds in the adiabatic limit to the dimeriza- 
tion in a Heisenberg model with alternating couplings. 
Its average over all sites I, 



1 N 



N 



(6) 



represents the uniform phonon displacement ,-and is 
closely related to the effective spin-spin couplingEm 



Jeff = J(l + 5). 



(7) 



It is well-known that in spin-Peierls substances like 
CuGeC>3 the transition from the uniform to the dimerized 
phase occurs at finite temperatures, here however, we 
consider a strictly one-dimensional model with a dimer- 
ization transition at T = 0. So to encounter dimerization 
one needs to go to temperatures as low as possible. 

By studying <5; at low temperatures in the parameter- 
range < ft/ J < 3 and < g < 2 we can clearly distin- 
guish two different phases of the model. We find cither 
pronounced periodic patterns in 5i (see Fig. ^|) or ran- 
dom fluctuations with no perceivable periodicity down to 
temperature T/J — 0.05. We note that as T is risen to 
higher temperatures the periodic patterns weaken and at 
higher T domains of periodic fluctuation are perceivable, 
which reflect a "reminiscence" of the periodic fluctua- 
tions present at T — (see second and third graph of 
Fig. |. 

It is of interest to see how the periodic patterns in the 
phonon displacement are connected with the formation 
of spin-dimers. So to complete the picture we calculated 
the local spin-spin correlations (ai ■ t?/+i). We find that 
the periodic patterns in the local displacement Si are well 
in phase with the formation of spin singlet states, which 
demonstrates the coupling of the phononic degrees of 
freedom and the spin degrees of freedom. This is also 
an a posteriori check of the consistency of our combined 
spin and phonon updates, (see the lowest three graphs of 
Fig. ||. As can be seen from the Fig. || the dimerization 
is not complete, since the correlation between the singlet 
states is not zero.) 

Thus Si proves to be a valuable tool, which allows more 
easily than the susceptibility to decide, for fixed param- 
eter values whether the model is in a dimerized state. 

We use this property to determine the approximate 
phase boundaries in the g — Q plain. To be concrete, 
we calculated Si for a 64 site system at temperatures 
T/J = 0.4, . . . , 0.1 and inferred from the periodic fluctu- 
ations the phase diagram shown in Fig. M. Since for all 
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points Si either shows random patterns or the periodic 
fluctuations, we can unambiguously determine, whether 
the model at a certain point in parameter space is in 
a dimerized phase. We cannot however state definitely, 
that the model is gapless if we encounter no dimerization, 
because we are limited both in system size and temper- 
ature. So rigorously the separation line between the two 
phases in Fig. ^| gives the minimal extension of the dimer- 
ized phase, i.e. in the thermodynamic limit the dimerized 
phase will be larger and the separation line in Fig. ^jwill 
be shifted to the left. Comparing with the results of 
Fig. U we find that the point g = 1.5 and fl = J where 
from the susceptibility we found that the model must be 
in a massive phase shows clear patterns of dimerization, 
and the point g = 0.8 and Q = J, for which the finite 
size argument was not conclusive, turns out to be in the 
dimerized phase. 




FIG. 9. Phase-diagram of the model in the g-Sl plain. 
Circles mark points with dimerization, squares points where 
no dimerization was found for N — 64 and T > 0.1J. 
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B. Effective coupling J e g 

We turn now to the averaged quantity S and study the 
effective coupling J e g as given by Eq. (0). It is worth 
noting that both in the dimerized and the undimerized 
phase 5 is nonzero i.e. the spin-spin coupling is always 
shifted by the presence of the spin-phonon coupling. Our 
analysis shows that for all choices of parameters J e e is a 
slowly decreasing function of temperature (Fig. ^), this 
behavior has been noted before for g = 1/ (4\/2) and 
n = 0.2 J in Rcf. | 
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FIG. 8. Upper three graphs: local displacement 

(bj + h) for T = 0.3J, 0.4J,0.5J. Lower three graphs: 
spin-spin-correlation (07 • ffi+i) for the same temperatures. 
All graphs are for N = 64, M = 80, Q = 2 J, g = 1.5, and 
averaged over 10,000 spin updates. 



FIG. 10. Effective spin-coupling J e ff/J versus temperature. 

/From the phcnomenological point it is interesting that 
a lowering of J with temperature is in agreement with 
Ref. |2l| where a comparison of inelastic neutron scatter- 
ing data on CuGeC>3 with calculated dynamic structure 
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factors S(q, f2) gave evidence for a decrease of J with in- 
creasing temperature. A similar behavior of J has been 
reported in Ref. p2[ 

With our results for J e ff we are in the position, to 
study the question to what extent J e g is an effective 
coupling in the proper sense, in other words, whether 
it can be used to rescale the thermodynamic quantities 
in such a way that the results for the Heisenberg model 
are regained. We find that the rescaled susceptibilities 
X(T) J e s(T) plotted versus temperature as shown in Fig. 
pi] do not fall together with the curves of the Heisenberg 
model. So one is forced to conclude that the effective cou- 
pling cannot be used to bring the curves with different 
parameter together. One cannot infer from this however, 
that essential new physics - maybe a sign of a frustration 
- is perceivable. The above result means only that J e g is 
not the correct variable to rescale the thermodynamical 
observables. 
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This scaling of the susceptibilities is of phenomenolog- 
ical importance also for CuGeC>3 , since it is-wpU known 
that the susceptibility of a frustrated modeJld'ta cannot 
be mapped by a rescaling of J onto that of the Heisen- 
berg model. T his len ds further evidence to the argument 
given in Sec. [II A that no agreement with the exper- 
imental susceptibility of CuGeC>3 can be achieved. By 
checking with data from the model with frustration we 
estimate that a frustration parameter of w 0.12 might go 
undetected in our analysis, because the statistical fluctu- 
ations hide the small effect on x f° r frustrations less than 
0.12. 
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FIG. 12. Rescaled susceptibilities x(J> 5> T, Q.) J^g/N ver- 
sus T/J o 9 ff 




FIG. 11. Rescaled susceptibilities x(J, 9,T,0,)J eB (T)/N 
versus T/J cff (T) 

On the other hand we find that to high precision the 
susceptibilities (for not too low T) can be rescaled by 



a global and temperature independent Jf ff /J (Fig. 12) 
That is we find 



x(J,g,T,n)j- 



off 



J X AFH (T/Jf ff ,Jf ff ) 



(8) 



where Jf ff is a global effective coupling , which we deter- 
mined by comparing the temperatures T max where the 
susceptibilities of the Heisenberg model and the model 
Eq. (|l|) have their maxima 



T 

jg _ r ± max 

"'eff — J 



yAFH 



(9) 



. £1 = 2J; T = 0.3J 
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0.0 



• g = 1.0; T = J 



0.0 0.4 0.8 1.2 1.6 2.( 

J/£2 



We note that there is a discrepancy to Ref. || where an 
overall (/-Factor was required to obtain a similar scaling 
law. 



FIG. 13. Displacement <5 = J c s/J— 1 versus g 2 (left side) 
and versus J/Q (right side) for N — 64 and Trotter number 
M = 80. 
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As a last point we show that S is proportional to g 2 /ft 
in the temperature range T < 2.0J, (see Fig. [l3| The 
constants of proportionality for the data shown in Fig. [l3| 
are c(T = 0.3J) = (2.38 ± 0.04)J and c(T = J) = 
(2.28 ± 0.02) J.) At higher temperatures we found de- 
viations from the simple law. This result is in qua. 
tative agreement with observations by Wellein et al 
who found a lowering of the dimerization with grow- 
ing frequency at 21 = 0, it also confirms calculations of 
Pytteo and UhrigEj, who found the spin-phonon Hamil- 
tonian (for translationally invariant phonons) equivalent 
to an effective spin Hamiltonian with spin-spin coupling 
oc g 2 /Q. 



C. Correlations 



With our method it is straightforward to calculate 
spin-spin-correlation functions < Si ■ Sj > as a function 
of temperature. We postpone a fit of the correlation func- 
tions in the strongly dimcrizcd phase and a detailed dis- 
cussion of the critical exponents in the gapless phase to 
a later publication and highlight at the end of this work 
the most spectacular features of the correlation functions 
only. 



10" 



10"' 



10" 



10 



oa = o.5J 

• £2 = 1.0J 

□ n = 2.0J 

• n = 4.oj 

a £2 = 20.0J 



8 5 ft 



8 12 
distance |i-j| of two sites 



16 



20 



FIG. 14. < Si ■ Sj > 
for g = 1.0 for a 64 site chain, 
fi = 0.5J, . . . , 20J at T = 0.3J. 



, as a function of distance 
M=80, phonon frequency 



(i) We find that the spin-spin correlations as a function 
of distance show a pronounced steplike behavior (see Fig. 
|l4] ) in the gapped phase. This is again a signature of 
the dimerization of the ground state. (Note that for a 
completely dimerized state we would have one step only, 
however, we find a finite correlation between the dimers.) 

(ii) For increasing Q at fixed g one finds the correlations 
more and more unaffected by the spin phonon coupling, 



which is due to the decoupling of spins and phonons as 
we observed before when studying the susceptibility. 



10" 



10 



s 10 



10 



og = 0.0 
• g = 0.2 

□ g = 0.8 

■ g = l.o 

g = i.i 

»g = 1.2 



i 2 



y • 

* I s • 

* ■ □ » 

Mas, 

' 5 i § 



■ • 

A A ■ f 



4 8 12 16 

distance |i-j| of two sites 



J 
20 



FIG. 15. Spin-spin correlation < Si ■ Sj > (— 1)* 3 , as a 
function of distance for g — 0, . . . , 1.2 for a 64 site chain, 
M=80, phonon frequency Q = 2 J at T = 0.3J. 

(iii) As a function of g we observe a crossover between 
the high-temperature and the low-temperature range. At 
low temperatures the correlations decrease with growing 
spin-phonon coupling (Fig. |l5| ). This is analogous to the 
dimerized model where one finds a lowering of the correla- 
tion with growing dimerization. For high temperatures, 
however, the inverse behavior occurs: the correlations 
for g = 1.5 are highest and for g = 0.2 are lowest. (For 
fl = 2 J the transition between low and high-temperature 
region occurs at T = (0.55 ± 0.05) J.) 




HeisenbBrg model 
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■ 9-1.0 my 




■ B 1 ^ 
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Heisenberg model fl^ 
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Q - on= 4J 
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FIG. 16. Correlation length £ in lattice units a for Q = 2 J 
(upper figures) and g — 1.0 (lower figures). The chain length 
is N = 64 and the Trotter number M = 80. 
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(iv) In the gapless region of the model the correlations 
decay exponentially with a correlation length that coin- 
cides for g < 0.2 with the known Heisenberg result 

£(T) = cT- v with v = 1 for N -> oo (10) 

(see Fig , |l6|b). For larger g we encounter deviations from 
Eq. (|l(j) which can best be fitted by allowing v to differ 
from 1. (e. g. v = 0.92 for g = 0.5 and = 2 J). In 
Fig. [l6|a,c we show £ as a function of g and A. We find 
that for temperatures 0.5 < T/J < 1 the correlation 
length £ decreases slightly with g and A. However, a 
more careful investigation is needed to give conclusive 
evidence for the functional dependence of f2, T) = 

(v) The pronounced maxima of the static structure fac- 
tor S(q, T) at q = tt for g = as shown in Fig. [l7] orig- 
inate from the logarithmic divergence of S(q,T = 0) of 
the isotropic Heisenberg model. (For a study as to what 
extent this divergence can be seen at finite T we refer 
to Ref. H|.) For the model Eq. (|l|) in the gapless phase 
the peak in S(q) at q = n is only slightly diminished 
as compared to the Heisenberg model, indicating that 
the models show a similar divergent behavior. On the 
contrary in the gapped phase the structure factor shows 
a broad maximum hinting toward a regular behavior of 
S(q = n,T = 0). 




0.2 - A £Cl= 4J;g=1.0 

* * g=0.0 



2.4 2.6 2.8 3.0 3.2 

q 



FIG. 17. Structure factor S(q) versus momentum q. 

V. CONCLUSIONS 

We have developed a modified quantum Monte Carlo 
loop algorithm to study the finite-temperature proper- 
ties of the one-dimensional isotropic antiferromagnetic 
Heisenberg model coupled to Einstein phonons. We have 
investigated in detail the magnetic susceptibility, specific 



heat, phonon occupation number, dimerization, and spin- 
correlation for chains with 16 up to 128 sites down to 
temperatures T = 0.05 J. 

We found evidence for a phase transition from a gap- 
less to a massive phase and used the periodic fluctuations 
in the local displacement to determine the approximative 
phase diagram of the model. Furthermore, the dimeriza- 
tion patterns we encountered were accompanied by the 
formation of spin singlet states. 

The temperature-dependence of the effective spin-spin 
coupling J e g was calculated and it was shown that 
Jc ff (T) = (1 + c(T)g 2 /il)J, for T/J < 2.0. 

At low temperatures the mean phonon occupation 
number was found to be proportional to the square of 
the spin-phonon coupling constant and inversely propor- 
tional to the square of the phonon frequency. 

We demonstrated that the susceptibility of CuGeC>3 
cannot be described by the model. 
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VII. APPENDIX 

The Trotter decomposition of the partition function Z 
we used as basis of our quantum Monte Carlo method 
reads: 

Af-l 

Z= Jim V TT 

fco,—,fc3M-l 3=0 

ffl <*W,i|e- /W(ij6,)/J, 'l«!w+i.« >) 

JV/2-1 

( II <h j+ i,2i\e-^ /M \k 3j+2 , 2l > 

1=0 

< fc3 J+ 2,2/ + l|e- /3H2 ' + l/M |fc3, + 3,2; + l >)■ (11) 

k denotes the spin and phonon configurations of the lower 
and upper edge of the interacting plaquettes, the inverse 
temperature is f3, the Trotter number M, and the Hamil- 
tonian for the spin-spin and the spin-phonon interaction 
is 

Wi = ^(^ +1 -l)(l + +&,)). (12) 







1 E. Pytte, Phys. Rev. B 10, 4637 (1974). 

2 M. Hase, I. Terasaki, and K. Uchinokura, Phys. Rev. Lett. 
70, 3651 (1993). 

3 D. Augier, D. Poilblanc, E. Sorensen, and I. Affleck, Phys. 
Rev. B 58, 9110 (1998). 

4 D. Augier, P. Hansen, D. P oilblanc, J. Riera, and 



E. Sorensen, |cond-mat/9807265 . 

5 G. Wellein, H. Fehske, and A. P. Kampf, Phys. Rev. Lett. 
81, 3956 (1998). 

6 A. Weifie, G. Wellein, and H. Fehske, |cond-mat/9901262 . 

7 D. Augier and D. Poilblanc, Eur. Ph ys. J. B 1, 19 (1998 ). 

8 A. W. Sandvik and D. K. Campbell, |cond-mat/990223C . 

9 A. W. Sandvik, R. R. P. Singh, and D. K. Campbell, Phys. 
Rev. B 56, 14510 (1997). 

10 R. Werner, C. Gros, and M. Braden, |cond-mat/9810038 . 

11 H. G. Evertz, G. Lana, and M. Marcu, Phys. Rev. Lett. 
70, 875 (1993). 

12 H. G. Evertz, in Numerical Methods for Lattice Quantum 
Many-Body Problems, edited by D. J. Scalapino (Addiso n 
Wesley, Reading), to be published (= cond-mat/9707221 ). 

13 G. Castilla, S. Chakravarty, and V. J. Emery, Phys. Rev. 
Lett. 75, 1823 (1995). 

14 J. Riera and A. Dobry, Phys. Rev. B 51, 16098 (1995). 

15 K. Fabricius, A. Khimper, U. Low, B. Biichner, G. 
Dhalenne, T. Lorenz, G. Dhalenne, and A. Revcolevschi, 
Phys. Rev. B 57, 1102 (1998). 

16 A. Khimper, Z. Phys. B 91, 507 (1993). 



17 R. J. Bursill, R. H. McKenzie, and C. J. Hamer, |cond- 
mat/9812409] 



^R. E. Rapp, E. P. de Souza, H. Godfrin, and R. Calvo, J. 
Phys.: Condens. Matter 7, 9595 (1995). 

19 M. Weiden, J. Kohler, G. Sparn, M. Koppen, M. Lang, 
C. Geibel, and F. Steglich, Z. Phys. B 98, 167 (1995). 

20 R. Fehrenbacher, Phys. Rev. B 49, 12230 (1994). 

21 K. Fabricius and U. Low, Phys. Rev. B 57, 13371 (1998). 

22 G. S. Uhrig, Phys. Rev. Lett. 79, 163 (1997). 

23 G. S. Uhrig, Phys. Rev. B. 57, R14004 (1998). 

24 K. Fabricius, U. Low, and K.-H. Mutter, Phys. Rev. B. 51, 
8270 (1995). 



10 



